Entanglement oscillations in open Heisenberg chains 
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We study pairwise entanglements in spin-half and spin-one Heisenberg chains with an open bound- 
ary condition, respectively. We find out that the ground-state and the first-excited-state entangle- 
ments are equal for the three-site spin-one chain. When the number of sites L > 3, the concurrences 
and negativities display oscillatory behaviors, and the oscillations of the ground-state and the first- 
excited-state entanglements are out of phase or in phase. 
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I. INTRODUCTION 

Recently, the study of entanglement properties in 
Heisenberg systems has received much attention 0- [28$ . 
Quantum entanglement is one of the most intriguing 
properties of quantum physics |29i l3Cj and the key ingre- 
dient of the emerging field of quantum information theory 
and processing 31], and can be exploited to accomplish 
some physical tasks such as quantum teleportation [32l ]. 

In most of the previous studies on entanglement of 
many-body states, a periodic boundary condition is as- 
sumed for spin chains. On the other hand, spin chains 
with an open boundary condition (OBC) have been used 
to construct spin cluster qubits [33T.l34| for quantum com- 
putation and employed for quantum communication from 
one end to another [3^, |3(| . Perfect state transfer has 
been obtained via the open chain without requiring qubit 
coupling to be switched on and off |3o| . These investi- 
gations reveal that open chains are of great advantage 
in implementing quantum information tasks. Thus the 
study of the entanglement structure in open spin chains 
will be of importance as the entanglement underlies oper- 
ations of quantum computing and quantum information 
processing. 

Most of the systems considered in precious studies are 
spin-half systems as there exists a good measure of entan- 
glement of two spin- halves, the concurrence |37j . which 
is applicable to an arbitrary state of two spin halves. On 
the other hand, the entanglements in higher spin sys- 
tems are not well-studied due to the lack of good opera- 
tional enta ngle ment measures. Here, we will use the neg- 
ativity |38l l39| | to investigate entanglement in spin-one 
systems. In this paper, by using the concept of concur- 
rence or negativity, we study pairwise entanglement in 
spin-half and spin-one Heisenberg chains with an OBC, 
respectively. 

The paper is organized as follows. In Sec. II, firstly, we 
give exact results of the ground-state and first-excited- 
state pairwise entanglements for the three-qubit and 
four-qubit spin-half open Heisenberg model; and then, 
we present numerical results of the corresponding entan- 
glements for the 5 ~ 10-qubit open Heisenberg model; 
finally, we investigate open boundary effects on thermal- 
state pairwise entanglement for 2 ~ 6-qubit Heisenberg 



model. In Sec. Ill, we study pairwise entanglement in 
spin-one open Heisenberg chains. We conclude in Sec. IV. 



II. SPIN-HALF SYSTEM 

The Heisenberg Hamiltonian for the chain of L qubits 
with an OBC is given by 

L-1 

H = J ( 2 ^ ■ Si+i + g) 

i=i 

L-1 

= 5>S M+ i, (1) 

1=1 

where Si is the spin-half operator for qubit i, iSj.i+i = 
4(1 + Si ■ <Xj+i) is the swap operator between qubit i 
and i + 1, and ai — (o~i x , <?iy, <?iy) is the vector of Pauli 
matrices. In the following discussion, we assume J = 
l(antiferromagnetic case). 

Due to the SU(2) symmetry in our Hamiltonian, the 
concurrence quantifying the entanglement of two qubits 
is given by [4(J, |41j 

dj = max{0, -2(S< ■ Sj) - 1/2} = max{0, -<<%)}, (2) 

we see that the entanglement is determined by the ex- 
pectation value of the swap operator. 

In the three-qubit case, the Hamiltonian can be written 

as 

H = S\2 + 1S23 

= 2(Si • S 2 + S 2 • S 3 ) + 1 

= (Si + s 2 + S3) 2 - s 2 - s 2 - s 2 - 2S a • s 3 + 1 

= (S 1 + S 2 + S 3 ) 2 -(S 1 + S 3 ) 2 -S 2 + 1, (3) 

because (S1 + S2 + S3) 2 , (S1 + S3) 2 and S 2 . commute with 
each other, we can use the standard angular momentum 
coupling theory to calculate all the eigenvalues of this 
system: firstly, Si couples with S3, then they couple with 
S 2 again. The results are 

E = -l(2), E 1 -1(2), E 3 = 2(4), (4) 
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where the number in the bracket denotes the degeneracy. 
The ground state is analytically given by 



L=3 



L=4 



= ^=(|ooi)-2|oro) + |roo)), 



or I* 



<2)\ 



1 

i 



(|110> - 2|101> + |011», (5) 



and the first-excited state is 



i(|001) - |100)), 

"- = ^(ino) -ion)). 



(6) 



Because the three-qubit Hamiltonian has an exchange 
symmetry, namely, the Hamiltonian is invariant after 
swapping qubits 1 and 3, (S12) = (S23) = \ (H) = -j. 
Thus, from Eqs.© and Q, the concurrence of two qubits 
in the ground state is found to be 

C° 2 = 1/2, (7) 

and for the first-excited state, the concurrence is 

C\ 2 = 0. (8) 

We see that the ground state is an entangled state, 
whereas the first-excited state is not. 

Now we consider the four-qubit case, the Hamiltonian 

is 



H — S12 + S23 + S34, 



(9) 



from which, we obtain all the eigenvalues of this system 
as follows 

E = -V3(l), £i = l-V2(3), #3 = 1(3), 
£ 4 = V3(1), £ 5 = 1 + V2(3), #6 = 3(5), (10) 
Then, the ground state is 
1 



l*o) 



|0011) - (2 + V3)|0101) 



V24+12V3 
+ (1 + \/3)|0110) + (1 + V3)|1001) 
-(2 + V3)|1010) + |1100)], (11) 



and the first-excited state is 



or |tf< 2) ) 



1 



— 10001) + (1 + %/2)|0010) 



V8 + 472 
-(1 + \/2)|0100) + |1000)] 



= [|0011)- (1 + v / 2)]0101> 



V8 + 4V2 
+ (1 + V2)|1010)- |1100)], 



or \V\ 



(3)\ 



-|0111> + (1 + \/2)|1011) 

(12) 



V8 + 472 

-(l + \/2)|1101) + |1110)] 




FIG. 1: The ground-state (star line) and first-excited-state 
(square line) nearest-neighbor concurrences versus site num- 
ber for L — 3 ~ 10 in the spin-half open Heisenberg chains. 



Thus, from Eqs.©, ifTT^l and (fT2"|) . we get the concur- 
rences of the ground state as 

C° 2 = C7 3 ° 4 = 3 + 2V l = 0.8660, C7 2 ° 3 = 0, (13) 
3 4 + 2^3 

and for the first-excited state, the concurrences are 

0^=01=0, ^3 = 1^ = 0.7071, (14) 

From the analytical results of the concurrences of the 
ground state and the first excited state, we observe that 
the concurrence oscillations emerge. 

Then we calculate the nearest-neighbor concurrences 
for the cases L = 5 ~ 10 numerically, and the results are 
shown in Fig. 1. From Fig.l, we see that the ground-state 
and the first-excited-state concurrences oscillate when 
the site index increases, and they are out of phase with 
each other when L > 3 (in the case of L = 10 , the up- 
ward convex of the square line in the region of i — 1 ~ 3 
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FIG. 2: The thermal concurrences between qubits 1 and 2 in 
the 2 ~ 6-qubit spin-half open Heisenberg chains. 



and 7 ~ 9 is not so obvious). The reason of oscillatory 
behaviors is as follows: for spin 2, there is a competi- 
tion between spins 1 and 3, and they both favor being 
maximally entangled with spin 2. If spin 2 shares a large 
entanglement with spin 1, then it is less entangled with 
spin 3, and vice versa. Thus the oscillatory feature ap- 
pears. 

Now, we study entanglement in the thermal state, 
namely, consider the finite-temperature case. The state 
of a system at thermal equilibrium is described by the 
density operator p(T) — exp(—f3H)/Z, where (3 — 
X/ksT, ks is the Boltzmann's constant, which is assumed 
to be 1 throughout the paper, and Z = Tr{exp(— (3H)} 
is the partition function. The entanglement in the ther- 
mal state is referred as thermal entanglement. Due to 
(Sij) = Tr[6>ij • p], then from Eq.Q, the thermal concur- 
rence quantifying the thermal entanglement is given by 



Cij (T) = max{0, -Tt[5 y • p(T)}}, 



(15) 



by using the above equation, the numerical results of the 
thermal concurrences for L — 2 ~ 6 are shown in Fig. 2. 
We can see that the thermal entanglements for even 
qubits are larger than those for odd qubits. At a fixed 
lower temperature, for even L the thermal entanglements 
decrease as the number of qubits increases, whereas for 
odd L, on the contrary, the entanglements increase as 
L increases. As L is very large, the effects of parity of 
L vanish. The threshold temperatures T t hS are indepen- 
dent of the qubits' number when L is large, and when the 
qubits' number is toward to infinity, they will approach a 
point, which is estimated as Tth = 1-71658 ~ 1.71659i"T. 



III. SPIN-ONE SYSTEM 

The Heisenberg Hamiltonian for the chain of L spins 
in spin-one system of an OBC is 



L-l 



H 



n+x, 



(16) 



i=l 



where Si is the spin-one operator for spin i, J is the 
exchange constant. Here we also set J = 1. 

For the case of higher spins, a non-entangled state has 
necessarily a positive partial transpose (PPT) according 
to the Peres- Horodecki criterion 38] . In the case of two 
spin halves, and the case of (1/2,1) mixed spins, a PPT 
is also sufficient. However, in the two spin-one particles, 
a PPT is not sufficient in general. The two-spin state 
here displays a SU(2) symmetry, and for such a state, 
the PPT condition is necessary and sufficient 0> C3- 
The quantitative version of the criterion was developed 
by Vidal and Werner j3^. They presented a measure of 
entanglement called negativity that can be computed ef- 
ficiently, and the negativity does not increase under local 
manipulations of the system. So, we may use negativity 
to exactly characterize the two-spin entanglement prop- 
erties of our system. 

For L = 3 case, we can write the Hamiltonian that 



H 



Si • S2 



- S2 • S3 

S 2 + S 3 ) 2 -(S 1 + S 3 ) 2 -S 2 .], (17) 



using the similar method in the above section, all the 
eigenvalues of this system are directly obtained as 



E = -3(3), E l = -2(1), 
£4 = 0(3), E 5 = 1(5), 

The ground state is obtained as 



E 3 
E 6 



-1(8), 
2(7). (18) 



1 



I* 



(2)\ 



(|002) + 6|020) + |200)- 3|011) 
60 

2|101> - 3|110», 

' (— 2|012> + 3|102) +3|120> + 3|021) 



I* 



(3)\ 



60 



+2jl21) - 3|211>), 

and the first-excited state is 
1 

+ 1021 



3|201> — 2|210) — 4|111», 
' (|022) + 6|202) + |220) - 3|112) 



(19) 

(-|012) + 1 102) - 1 120) 

|201) + |210)). (20) 



In spin-one system, the negativity quantifying the pair- 
wise entanglement is 



Mi 



1 



max[0, (Sij) — 
- max[0, — (S^ 



(S 7 ; • So 



(21) 



FIG. 3: The ground-state (star line) and first-excited-state 
(square line) nearest-neighbor negativities versus site number 
for L = 3 ~ 6 in the spin-one open Heisenberg chains. 



FIG. 4: The thermal negativities between spins 1 and 2 in the 
2 ~ 6-spin spin-one open Heisenberg chains. 



where £y = S, • Sj + (Si • Sj ■ ) 2 — I is the swap operator be- 
tween the spin-one particles i and j . For the Hamiltonian 
of L = 3, the spins 1 and 3 have an exchange symmetry 
which leads to (Si • S2) = (S2 • S3) and (S12) = (S23}, so 
we have A/12 = A/23. 

Thus, from Eqs.(19) and (20), the expectation values 
of S12 and Si • S2 in the ground state are given by 

(S U )° = ± (Si.S 3 -)° = -| (22) 

and in the first-excited state, they are 

{S ia ) 1 = -1, {Si ■ Sj) 1 = -1. (23) 

Then, from Eq.(21), we obtain the corresponding nega- 
tivities as 

A/? 2 = l/3, ^ = 1/3. (24) 

From the above equation, we can see that the negativities 
in the ground state and the first-excited state are equal. 
When L > 3, the negativities in the ground state and 
the first-excited state are not equal, so this equality is a 
mathematical accident. In addition, because the reduced 
density matrices in the ground state and the first-excited 
state arc different. 

Numerically we calculate the negativities of L = 4 ~ 6, 
the results are given in Fig. 3. We can see that the neg- 
ativities oscillate when L > 3, and the ground-state and 
the first-excited-state negativities are out of phase with 
each other for L — 4 and 5 cases, but they are in phase 
for L = 6. Then from Fig.f and Fig. 3, considering the 
problem of the oscillatory phase of the ground-state and 
the first-excited-state entanglements, we find out that: 
on the edge of the Heisenberg chains, due to the effect of 
the boundary condition, the ground-state entanglements 



are biggest, thus the oscillatory phases of the ground- 
state entanglements are independent of the spins' num- 
ber and invariant; but for the first-excited-state entangle- 
ments, the entanglements on the edge are small when the 
spins' number is not big, then because of the existence of 
the competition, the nearest entanglements will be large, 
thus the ground-state and the first-excited-state entan- 
glements are out of phase with each other; along with the 
spins' number increasing, the first-excited-state entangle- 
ments on the edge will be larger in contrast with the cor- 
responding ground-state entanglements, which leads to 
the nearest first-excited-state entanglements diminishing, 
thus when the chain has enough spins, the ground-state 
and the first-excited-state entanglements are in phase. 

Becase (S, • Sj) = Tr[(S, • S,-) ■ p],((S< ■ S,) 2 ) = 
Tr[(Sj ■ Sj) 2 ■ p], then from Eq.(2I), the thermal nega- 
tivity quantifying the thermal entanglement is obtained 
as 

A/-,(T) = \ maxj^Ti-KS, ■ S,) 2 • p(T)] - 2} 

+i max{0, 1 - Tr[(Si • S,-) • p(T)] 
-Tr[(S, • S,) 2 • p(T)]}, (25) 

according to Ea.l|25[). we calculate the thermal negativi- 
ties of L — 2 ~ 6 numerically, which are shown in Fig. 4. 
Comparing Fig. 4 with Fig. 2, we find out that the behav- 
iors of thermal entanglement in spin-one system is similar 
to those in spin-half system, but the threshold temper- 
atures in spin-one system are all smaller than those in 
spin-half system, and when the number of spins tends to 
infinity, they will approach a point which is in the region 
of T th = 1.26753 - 1.26758JT. 
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IV. CONCLUSION 

In the above, we have studied pairwise entanglements 
in spin-half and spin-one Heisenberg chains with an OBC, 
respectively. Some analytical and numerical results of 
the ground-state and first-excited-state pairwise entan- 
glements in the chains for several spin particles have been 
presented, and we find out some interesting results: the 
ground-state and the first-excited-state negativities are 
equal when L = 3; and when L > 3, the concurrences and 
negativities both oscillate, which results from the compe- 
tition between the two spins on both sides of each spin. 
We also have given numerical results of the thermal-state 
pairwise entanglements in the chains for a few spins, the 
results reveal that the thermal entanglements in the two 
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